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(5) INTRODUCTION 


The goal of this research is to improve ultrasonic classification of breast lesions and guide 
decisions regarding biopsy requirements, especially for small lesions and those in young, dense 
breasts, which are particularly difficult to evaluate with mammography. The research is designed to 
go beyond current capabilities by developing a set of complementary ultrasonic morphological 
analysis procedures (UMAPs) that analyze digital ultrasonic echo data. Each UMAP is optimized 
for extracting a particular quantitative lesion feature that describes the spatial configuration or 
acoustic properties of the lesion. The set of complementary UMAP features will be analyzed with 
statistical procedures to derive reliable and objective lesion classification for use in patient 
management. 
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(6) BODY OF ANNUAL SUMMARY 


The research conducted during the first year of this program addressed issues itemized in the 
Statement of Work of our proposal. The research is planned to derive quantitative UMAP 
evaluations of features that are now subjectively described; these include “echogenicity,” 
“heterogeneity,” “shadowing,” and lesion boundary characteristics. Progress during this report period 
will be described in terms of the constituent steps needed to evaluate our concepts. These include: 
ultrasonic parameter ima^ng; region-of-interest definition; implementation of complementary UMAP 
techniques; and, initial classification efforts. 

During this report period, we developed quantitative ultrasonic parameter imaging techniques 
to process original digital data and display cross-sectional images of B-mode ultrasonograms and 
selected ultrasonic spectral parameters. The clinical data being processed was obtained previously 
at three clinical sites using ATL Ultramark 9 systems with linear LI0-5 arrays. Data sets include 
digitized radio-frequency (RF) echo data acquired at the transducer, before any conventional 
processing was applied. Calibration signals from each system are also stored for use in removing 
system-dependent signal artifacts. Complete clinical records (including subsequent biopsy reports) 
are stored anonymously for each of 40 biopsy-proven cancerous lesions and more than 100 biopsy- 
proven non-cancerous lesions. Also stored are the clinicians’ level-of-suspicion that each lesion was 
cancerous; these levels were based on standard interpretation schemes applied to conventional 
ultrasonograms. Digitized RF data records were obtained for several complete scans of each lesion 
as well as ultrasonic scans of the contralateral breast. 

The new digital algorithms we developed during this report period use MATLAB applications 
for generating cross-sectional tissue images. For B-mode imaging, stored time-gain-control (TGC) 
data are used to compensate for operator settings and overall system gain. Stored RF data are 
processed with an analytic-signal-magnitude routine to obtain high-quality video signals, which are 
employed for B-mode images. 

During this period, we implemented and tested spectral processing techniques that involve a 
sequence of procedures to compute and display frequency-dependent features of tissue backscatter. 
After TGC compensation, RF echoes along each scan line are analyzed with a sliding 64-point Fourier 
Transform algorithm (incorporating a Hamming weighting function). At each location spectral 
amplitudes are computed in dB over a selected frequency range (5-9 MHz). To extract summary 
spectral features, linear regression is applied to compute the spectral intercept (dB; extrapolation to 
zero frequency) and midband fit (dB; value of the regression line at the center frequency). 

Local intercept and midband fit values are next calibrated with respect to an absolute 
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reference level, using stored calibration data obtained for each clinical system. Our processing 
includes two calibration procedures. The first removes frequency characteristics of the transducer 
and pulser-receiver; here, relevant data were obtained with planar focal-plane reflectors. The second 
procedure removes range and frequency effects of ultrasonic diffraction; relevant calibration data 
were obtained with rubber blocks embedded with diffuse distributions of 10-pm glass spheres. 

The final result of these procedures is a set of images and files containing corresponding 
values of; B-mode (video) signals; intercept values; and, midband fit values. 

During this period, the above procedures were all implemented and tested on clinical data 
using MATLAB. We also developed region-of-interest (ROI) procedures using algorithms that can 
be applied to any of the above image types. As described below, medical personnel must identify 
several types of ROI’s in each scan plane; thus, we developed simple user-interfaces to define, review, 
and edit manually selected ROI’s. The ROFs include: lesion boundaries; anterior segments (directly 
anterior to and lateral to the lesion); similar segments posterior to the lesion; lateral regions on either 
side of the lesion; and, the breast surface (for cases where stand-off coupling pads were used during 
scanning). 

After developing the above imaging and ROI algorithms, we implemented processing software 
for key UMAP functions. The following measurements were implemented: 

1) “Echogenicity”; This is now computed as the mean spectral intercept value within the 
outlined lesion boundary. Intercept was selected for this measure because it is not 
affected by typical attenuation losses in intervening tissues. 

2) “Heterogeneity”; This is now computed as the standard deviation of the spectral 
midband fit within the lesion. We also implemented several alternative measures of 
statistical dispersion for study during the remainder of this program. Midband fit 
values were chosen for this measure since they can be estimated with higher precision 
than other spectral descriptors. 

3) “Shadowing”: This parameter measures the attenuation in lesions relative to the 
attenuation in surrounding tissues. We have implemented software that quantitatively 
estimates lesion attenuation coefficients by computing the difierence between midband 
fit values in shadowed regions (directly posterior to the lesion) and values in similar 
unshadowed regions (laterally displaced from shadowed regions). The computation 
accounts for lesion thickness, determined from the traced lesion boundary. 
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To complement these features, we have now started to examine a number of procedures to 
quantitatively characterize lesion-boundary features such as lobulation and eccentricity. 

We have tested and refined the above spectral measurements on a selected subset of clinical 
cases. The current UMAP feature set, augmented by lesion-boundary features and comparisons with 
contralateral-breast features, will be applied to the entire set of clinical data during the next program- 
year. Classification will be quantified using discriminant analysis and Receiver Operator 
Characteristics (ROC). To determine incremental benefits for identifying cancerous lesions, we will 
compare ROC curves for UMAP classification with ROC curves for clinicians’ level-of-suspicion 
(based on conventional ultrasonic interpretation schemes). 
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(7) SUMMARY APPENDICES 


1) List of key research accomplishments 

• implementation of software for computing, calibrating, and displaying ultrasonic 
spectral features of breast tissues 

• implementation of software and simple user-interface for defining breast lesions 
and relevant regions-of-interest 

• implementation and testing of software for lesion feature extraction including 
quantitative measures of 

- echogenicity 

- heterogeneity 

- attenuation 

• initial specification of techniques to quantify the spatial configurations of lesions 

• initial testing of concepts on selected cases preparatory to comprehensive 
application to all cases and statistical evaluation of lesion identification, to be 
performed during the following period. 

2) List of reportable outcomes 

“Ultrasonic spectrum analysis procedures for breast cancer classification,” S.K. Alam, 

F.L. Lizzi, E.J. Feleppa, T. Liu and A. Kalisz, 24th International Symposium on 

Ultrasonic Imaging and Tissue Characterization, Arlington, VA, June 2-4, 1998. 

Several additional abstracts and manuscripts are in preparation or review. 

3) Copy of abstract (on following page) 
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Ultrasonic spectrum analysis procedures for breast cancer classification 

S.K. Alam, F.L. Lizzi, E. J. Feleppa, T. Liu and A. Kalisz 
Riverside Research Institute, 330 West 42nd Street, New York, NY 10036 


We have developed a series of spectrum analysis procedures designed to quantify ultrasonic breast 
cancer evaluations. The procedures have been planned to improve upon B-mode differentiation of benign 
and malignant lesions, which employs features such as “echogenicity,” “heterogeneity,” and “shadowing.” 
Our goal is to replace each of these subjective features with corresponding quantitative parameters based on 
spectrum analysis of radio-frequency (RF) echo signals in order to remove operator dependence and to 
permit objective discrimination. 

Our technique involves an image-based approach to calibrated spectrum analysis. This has been 
implemented using RF data digitally acquired from several clinical sites using an ATL Ultramark 9 system. 
The first step is the digital synthesis of spectral-parameter images derived with sliding-window Fourier 
transform techniques, as described in previous reports for other organs. Images of uncalibrated local values 
of spectral intercept and midband fit are generated using a new MATLAB® (The MathWorks, Inc., Natick, 
MA) implementation. These quantitative images are then calibrated using the spectrum of a planar target 
together with a range-dependent diffraction correction for each parameter. Diffraction correction employs 
power spectra measured from diffuse scatterers in a gel or rubber block; it depends upon the specific 
transducer array and transmit focal-length used in each examination. 

Classification parameters are derived after tracing the boundary of breast lesions on midband fit 
images; each parameter replaces a specific B-mode descriptor. “Echogenicity” is measured as the mean 
spectral intercept within the lesion, since this value is not significantly affected by frequency-dependent 
attenuation in intervening media. “Heterogeneity” is measured as the statistical dispersion of midband fit 
values within the lesion. This definition is motivated by previous analysis that demonstrated how the 
histogram of midband fit, and its variance, can be related to tissue homogeneity. “Shadowing” is quantified 
by measuring mean midband fit values in two comparable regions of posterior tissues that are shadowed 
and not shadowed, respectively, by the lesion. The lesion attenuation coefficient is estimated from the 
difference between these mean values and the lesion thickness. 

Initial results on biopsy-proven cases are promising. We are now investigating additional 
descriptors for lesion surfaces to quantify the “smoothness,” “lobulation,” and “invasiveness” categories 
that have proven useful in B-mode evaluations. 
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